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ABSTRACT 


Title  of  Tliesis: 

UMDFET  —  A  TWO-DIMENSIONAL  GENERAL  DEVICE  SIMULATOR 
AND  ITS  APPLICATION  IN  EPROM  ANALYSIS 

Name  of  degree  candidate:  Zezhong  Peng 
Degree  and  Year:  Master  of  Science  19S9 

Thesis  directed  by:  Jeffrey  Frey.  Professor  of  Electrical  Engineering  Department 


A  generalized  single-electron-gas  energy  transport  model  for  simulation  of  sub¬ 
micrometer  dimension  GaAs  devices  is  described,  in  which  the  effects  of  the  heat 
flux  and  the  potential  energy  of  the  electrons  in  the  upper  valley  are  included.  The 


electron  temperature  includes  the  effect  of  introduction  of  the  potential  energy. 
A  numerical  method  for  rigorous  solution  of  the  current  continuity  and  energy 
transport  equations  is  presented.  A  two-dimensional  general  device  simulator  — 
UMDFET  was  developed  using  some  of  the  principles  of  the  model  and  has  been 


used  to  simulate  EPROM  devices. 
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Chapter  1 

General  Information 


UMDFET  (University  of  Maryland  Field-Effect  Transistor  Simulator)  [1,2, 3, 4, 
5,0]  is  an  newly  develoi)ed,  user-oriented,  two-dimensional  numerical  device  sim¬ 
ulation  program,  which  may  be  used  for  MOSFETs,  LDD  devices,  MESFETs, 
llEMTs  and  planar  multi-layer  structures.  The  device  material  can  be  Si,  GaAs 
or  any  other  material  for  which  a  user  can  provide  characteristics.  Both  a  drift- 
diffusion  (DD)  model  and  an  energy  transport  (ET)  model  are  available.  In  the 
latter,  the  two-dimensional  energy  balance  equations  are  rigorously  solved,  includ¬ 
ing  heat  flux  and  diffusion  terms,  and  impact  ionization  £is  well  as  the  potential 
energy  of  electrons  in  satellite  valleys  in  compound  semiconductors.  UMDFET  is 
capable  of  simulating  the  performance  and  reliability  of  devices  with  deep  submi¬ 
cron  channel  lengths  and  under  high  bias  operation. 
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1.1  User’s  Guide 

UMDFET  is  «  fricndli’  simulation  program  wliich  can  interact  with  a  user,  and 
guide  the  user  to  a  successful  simulation. 

UMDFET  has  the  following  capabilities: 

1.  With  the  energy  transport  model,  UMDFET  can  accurately  simulate  de¬ 
vices  with  effective  channel  lengths  below  0.1  micron.  Sucli  short-channel 
effects  as  velocity  overshoot  and  electron  temperature  lowering  arc  accu¬ 
rately  doscriljcd[  l,5].  Refer  to  TM  of  OPTION  card  in  Chapter  2. 

2.  UNIDFE']'  is  ahlo  to  calculate  hot  electron  induced  gate  injection  current 
and  substrate  current  [3,0].  Refer  to  IG  in  OUTPUT  card  for  gate  current; 
IMP  and  TM  in  OPTION  card  and  IMP  in  OUTPUT  card  for  substrate 
current. 

3.  User  can  ea.sily  specify  materials,  i.e,  Si  and  GaAs,  or  any  other  material 
whose  characteristics  arc  known.  Refer  to  LAYER  in  Chapter  2. 

4.  Capability  of  analysis  of  planar  multi-layer  device  structures  enables  user  to 
simulate  conventional  MOSFETs,  LDD  MOSFETs,  MESFETs  and  IIEMTs. 
Refer  to  DE\'1CE,  LAYER,  PROFILE,  IMPLANT  and  LDDIMP  cards  in 
Chapter  2. 

5.  User  can  get  on-line  information  on  the  UMDFET  simulation  process  by 
specifying  PRINTLEVEL,  so  he/she  can  check  whether  the  simulation  is 
processing  properly. 

6.  Convergence  at  very  high  bias  is  good.  Methods  of  global  convergence  check¬ 
ing,  2D  potential  extrapolation  and  convergence  refining  are  used  to  improve 
convergence[2,3].  Refer  to  OPTION  and  END  cards  in  Chapter  2. 


2 


7.  User  lias  good  conlrol  of  UMDFET  in  order  to  achieve  optimum  simulation 
on  a  specific  problem.  Refer  to  INTBIN,  NREFINE,  MCYCLE,  ERROR 
and  PRINT  in  Chapter  2. 

8.  UMDFET  has  two  kinds  of  simulation  bias  modes,  i.e.,  an  I-V  curve  gen¬ 
erating  bias  mode  and  an  EPROM  bias  mode.  Refer  to  BIAS  and  STEP. 

9.  UMDFET  can  use  previous  UMDFET  simulation  results  as  an  initial 
guess  for  a  now  simulation.  User  can  refine  the  convergence  of  a  previous 
simulation  to  any  desired  level,  so  that  a  hard-convergence  problem  can 
be  solved  incrementally.  Refer  to  INTBIN  and  NREFINE  in  OPTION  in 
Chapter  2. 

10.  User  has  three  options  for  mesh  generation: 

1)  user  defined  mesh; 

2)  taking  a  mesh  from  UMDFET.INT;  * 

3)  auto-mesh  generation. 

Refer  to  MESH  in  Chapter  2. 

11.  Auto-mesh  updating. 

12.  Versatile  doping  profile  generation. 

Refer. to  PROFILE  and  LDD  in  Chapter  2. 

13.  Good  interfaces.  Refer  to  UMDFET  File  System  and  Block  Diagram. 
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DEVICE  Simulation  for  O.I  micron  MOSFETS  -  Performance  and  Reliability 


Jettrey  Frey 

University  ot  Maryland,  College  Park.  MD  20742 


The  Energy  Transport  (£1 )  method  tor  solution  ot  the  Boltzmann  |  ransport  tquation  is 
an  efficient  and  economical  approach  to  the  simulation  of  the  performance  and  reliability 
of  devices  witn  active  lengths  down  to  only  hundreds  ot  Angstroms  1  nis  method  is  rapid 
enough,  even  with  accuracy-ennancing  measures  taken,  and  when  applied  to  two-dimensional 
problems,  to  be  used  at  an  maivioual  engineer  s  workstation 

While  original  ET  device  simulations  were  confined  to  use  ot  a  symmetrical  distribu¬ 
tion  function  current  simulations  using  this  method  either  aad  an  additional  moment  equa¬ 
tion,  to  account  for  neat  flow,  or  approximate  the  distribution  function  wiin  physically 
reasonable  functions.  We  have  adopted  the  latter  approach,  and  incorporated  tne  method 
into  our  two-dimensional  device  simulation  program  UMDFE  r,l  1  j  i  he  distribution  function  we 
use  Is  an  expansion  in  Legendre  polynomials  to  account  for  a  spnencaily  symmetrical  and 
an  asymmetrical  component,  ana  yields  energy  oistnoution  function  results  very  ciose  to 
those  obtained  by  MC  simuiaiions-at  a  very  small  fraction  ot  the  cost  in  addition  use 
ot  this  expansion  allows  the  E  l  moment  equations  fo  be  expressed  in  very  simple  and  easily 
calculable  form. 

We  have  used  these  methods  m  our  two-dimensional  E I  program  to  study  n-channel 
Silicon  MOSFET  s  with  gate  lengins  between  0  05  ana  0.25  microns  With  consiant-tieid 
scaling,  maximum  electron  temperature  in  these  devices  decreases  significantly  as  channel 
length  decreases  for  channel  lengins  less  than  0  t8  microns  This  pnenomenor..  which  is  a 
result  ot  the  decreased  ratio  of  electron  transit  time  to  energy  relaxation  time  as  gate 
lengths  are  shortened,  indicates  tn^t  snort-cnar'  el  devices  can  oe  more  reliable  tnan 
similar  long  channel  devices  as  indicated  also  by  recent  experiments  |2J 

Our  short-channel  ET  simulations  which  require  only  SUN-levei  workstations  show  that 
hot-electron  related  MOBFE I  reiiaoiiity  problems  should  become  less  threatening  for  silicon 
MOijFEl  s  with  very  short  channel  lengths  These  simulations  will  be  equally  usetui  tor 
predicting  basic  device  performance  factors  as  weii-if  suitaoie  pror -'s  reiaii  u  constants 
can  be  obtained--as  tor  predicting  nevice  reliability 


( 1  j  2  Peng  "UMDFET  A  Two-dimensional  General  Device  Simulator  and  ns  Application  in 
EPHOM  Analysis".  Master  s  Thesis  of  19B9.  in  the  University  of  Maryland 

(2J  G.G.  Bhahidi,  D  A  Antomadis,  and  h  I  Smitn.  "Reduction  ot  Cnannel  Hot-Eiectron- 
Generated  Substrate  Current  m  Sub- 1 50nm  Channel  Length  Si  MOSFE  T",  IEEE  Electron  Device 
Letters,  vol  9,  no.  10.  p  497.1988 


SUSCEPnBlUrY  of  culuum  ASSEMDE  AMPUFIEBS  to  SIN<JL£  pulses  of  intense 

MICROWAVE  Radiation 


Jobn  H.  hicAdoo,  W.  Mtdiari  BoDeii.  Robot  Seeley,  WiD  Caoe 
MLssdcu  Reseaicli  Corp^^ewiugton,  Va. 

Jeffiey 

Dcpamneiit  of  Ekctckal  Eiigbeoiiig 
Uozvetsity  of  Marylsid,  College  Pa±.  Md. 


Ganam  aisraiirie  (GaAs)  moaoMuc  microwave  mtegraied 
dicnats  (MMlCs)  are  staring  to  be  ioserted  into  fniHtary 
qrstoQs  heranse  they  ate  eiaiemely  compact.  light  wdght, 
aad  have  low  power  reqturr.menrs.  Berany  of  their  low 
noise  figure,  the  gnanasr  and  most  sensinve  of  these  de¬ 
vices  are  eacpccted  to  find  es^ially  wide  nse  as  rcceivcts 
for  radar  fieqoeocy  signals  and  herein  lies  a  problem.  The 
Sabmicroa  transistor  dements  have  low  beat  cqjadty  and 
dierefore  it  takes  little  energy  to  beat  them  to  point  of  per- 
manent  damage.  An  aiiciafl  employing  GaAs  MMIC  t^lar 
teceivexs  can  expect  to  have  the  leceiveis  damaged  by  a 
close  ap{Hoach  to  any  powerful  seatch  radar.  But  even 
woEcse,  the  incoiporadoo  of  these  devices  into  our  military 
systems  will  give  an  enemy  an  incentive  to  develop  elec¬ 
tromagnetic  directed  energy  systems  capable  of  destroying 

the  small  tTan<a'<atTr<;  frnfn  ran  of  leyml  Inlom^ifrg  FOT 

many  ^licatioos  it  will  be  pc^'bie  to  protect  cl^  receiv- 
ets  by  [facing  a  conventional  limiring  diode  between  the 
antenna  and  the  receivor.  However,  for  phased  anay  radar 
the  packing  den^  of  the  indtvidnal  receiver  elements  wiZl 
have  to  be  too  high  to  peonii  the  inclusion  of  macroscopic 
power  hmtting  diodes  and  therefexe  a  tnicroaxipic  solutiao 
fix  hardening  the  system  mna  be  foond. 

Hgore  1  shows  a  typical  experimeatal  seo^)  used 
fix  deienninmg  snscgeihiliiy  and  for  stndying  damage 
mechanisms.  The  circuit  is  used  to  measure  the  power 
xt^  to  the  device  nnder  test  (DUT),  die  ootpDt  power  and 
the  power  reflected  in  the  fundamental  ficqtgicy  and  several 
barmooics. 

Muiowave  radiatioa  fiom  a  (Erected  energy 
we^ion  is  expected  to  consist  of  a  ^gle  inieose  micro¬ 
wave  poise  lading  fiom  10  to  100  ns  at  1  to  12  GHz  with 
an  intensity  on  target  as  high  as  1  KwAan^  In  the  labora- 
kxy.  we  have  qmnlatftd  this  enviiomneat  by  iniecong  single 
poises  of  radjation  xnio  the  input  ports  of  nnpmteriRd 
devices.  The  testiug  procedme  to  determine  soscqttiblliiy 
was  smilar  for  all  devices.  First  a  low-ptTwor  pi^  was 
Bsed  to  voify  diat  the  device  was  fimc&ooing  oormaDy. 
Soccesstve  pidses  were  then  applied  wxdi  6dB  power  ineze- 
neacs  nndl  die  device  was  damaged.  Damage  was  asoer- 
tainrd  tfaroogb  varioos  tecfatQqoes.  In  the  c^y  tests,  the 
has  cmeots.  noise  figoie  and  small  signal  gain  were 
measned  after  each  poise.  A  cdonge  in  noise  figure  of 
dB  was  a  agoificant  iodicattg  of  damage  and  was  always 
aooompaoied  by  a  coae^tondingly  significant  redaction  in 
(Vnall  signal  g^an.  As  the  testing  proceeded,  it  became 
q:paieat  that  visual  inspectioQ  of  the  circuit  at  a  magnifi- 


Hgure  1 

A  typeal  experimental  seo^)  for  susc^idNliiy  tests  and 
hankoing  rslaied  studies. 

catkffl  of  about  SOCK  often  provided  an  even  more  strin¬ 
gent  test  for  damage  than  either  noise  fignre  or  qnall  sig¬ 
nal  gain,  and  theiefine,  in  the  later  tests,  we  nsed  only 
visml  rhfvJct  imd  changes  in  qgnat  pm  as  indica- 
tOES  of  damage. 

Fignie  2  ediibiis  the  main  resnlts  of  oor  snscep- 
EbOity  tests.  It  diows  the  device  type  vs  the  mcitleat  energy 
in  die  poise  diat  lesttUed  in  detectable  damage.  The  length 
of  (he  Mack  bats  is  propcxtioDal  to  the  average  of  aQ  the 


Incideat  Eaetgy  (nJ) 

Egnrc  2 

Tntjfiwit  eoef]gy  for  failtce  ftr  vadcxis  device  o^pes.  If 
(be  eneiOT  is  less  than  that  shown  by  the 

atriped  bars,  the  device  win  sorvive,  bot  if  it  is  azs  x 
than  the  bl^  bars,  (be  dence  win  oHtainly  fail  'Lhe 
are  io  oixkr  of  freoaeocy  match  between  the 
device  desigii  hegoeacy  and  w  test  ftegoeocy.  The 

Iasi  shftfg  for  the  coaesponding  device  type  tested,  while 
die  length  of  the  grey  bars  is  pttportioQ.'*!  to  the  aveiage  of 
an  (he  p»nnitin*ait>  diots.  If  ooe  applies  more  incident 
eoagy  fem  that  shown  by  &  dadc  bar.  dte  device  wQi 
ptofaably  fan.  bat  if  ooe  ipi:^es  les&  than  that  shown  by  the 
gity,  die  device  wfll  most  likely  survive.  'The  device  types 
arc  Atecviaifid  as  follows:  Powa=GaAsMESFET  power 
atspliCer,  LNA  s  low  noise  amplifier  employiDg  GaA$ 
aetal-semicondoctor  5eld  effect  transisinr  (MESFET, 
HEMT  K  low  nni^-  ymplifirr  employing  high  exction 
mobOity  uansistots.  Along  the  vctdcal  axis,  with  the  de¬ 
vice  types,  is  a  nember.  tj,  that  indicaies  the  fieqpency 
match  between  the  lest  freguency  and  the  design  fieqaeocy 
of  die  device, 

where/,  is  is  the  test  fioequcncy,/,  the  design  fieqneacy 
and  the  bandwidth.  The  devices  are  listed  in  Order  of 
beqaeticy  inaich.  with  the  worst  inatcfa  at  the  top.  Among 
the  low-noise  amplifiers,  one  can  discern  a  Uight  trend 
toward  greater  sascqxibility  as  the  fiequent^  match  im* 
pcovct.  The  lelanon  is  weak  because  (he  refiectivity  of  die 
toband  devices  increases  with  inoeasiDg  incident  power, 
wfaSe  that  of  the  om-of-band  devices  decreases.  T^isa 
iBEdz  of  the  power  dependence  of  traasiacr  mpedance. 
One  low  noise  aiqilifier  bodes  the  trend.  Tbisemploysa 
tranqanr  Aisigned  with  a  larger  immber  of  nar- 
xower  gales  than  its  Amedcan  coonieiparvs. 

Focossiag  CO  the  ovetall  picnire,  tbe  tests  diow 
^atGaAs  MMIC  and  KQC  devices  fail  at  stress  levels  that 
are  ooe  Older  of  magounde  less  Ibao  the  goals  set  by  miliary 


ptorm^  Psr«^”dy  wonisnme  ace  the  resale  oa  devices 
eo^'^ying  tngh  electron  naolaliiy  transigont  OuihfTs). 
These  are  mote  snscendble  by  an  order  of  magrritode  than 
the  <  iAr  r"**”l  field  cLTect  transistor  tecfaootogy. 

The  snscqtdbiii^  tpg?:  dtat  tnHifumft  ifac  problem 
were  largely  coo^eied  m  19S7.  Although  snscqxibQiQ' 
tests  we  coocmaing  as  new  devices  are  made  avashd^  oar 
effons  have  shifted  (owacds  research  aimed  at  hardening 
tcnlitfrmc  Diese  cniiirifws  fitQ  ixUD  two  caiegpncs:  tech- 
niqoes  that  the  power  trareamned  to  the  sensitive' 
MPjsbTn’  gstffx,  and  wdmiques  that  enable  the  gates  to 
wiftaant!  gicattr  abase.  Wc  have  been  fbensing  on  the 
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Damage  to  n  MESFET  gate  fiem  a  9Gaz.  lOns  pnlse 
at360W. 

bitifj-  category  by  soa^ring  damage  mechanisns.  Early 
qurfrcHf  pbotomiaogfqiis  of  damaged  gates  (Egnre 
3)  revealed  that  the  eoeegy  dgnsition  is  conoaurated  at 
gate  nrrtK  These  results  soggest  that  the  corneg  of  the 
aooice  and  drain  deettodes  dioald  be  xooaded  to  ledoce 
eiagrrir  field  ronqauraiioos  and  badafi  resistance  should 
be  added  to  spread  tbe  eoeigy  mcee  evenlbr  over  die  vadons 
pmVn  gpie  rhantu4c 

More  recent  wodc  eoplcTmg  specnosoopy  and 
high  ^peed  pbotogE^bs  of  developing  qaiks  has  revealed 
a  Mliire  mode  called  recdficanon  bahne  conunoD  to  GaAs 
MESFETs.  The  mecbaniaai  remits  in  the  destmetiaa  of 
the  rectifying  piopeny  of  the  Schotlky  joacdoos.  Ajnnc* 


boD  that  has  so  £aQt^  cao  do  loDger  biock  the  Qow  d 
cuLca  in  ritfaer  diiccdon  and  acts  as  as  ordioaiy  ohmic 
nxiflcr.  Catastrophic  iaHutc  follows  hem  coneat  ddven 
by  bias  volmges  through  the  bsiled  juncDon. 

Ftevioos  investigatois  employing  trains  of  lower 
power  poises  had  identified  metal  nrigmfioa  as  a  key  fitil- 
ore  mode  for  extended  eaposuic  to  modoate  levels  of 
miaowave  radiation’'^.  Qjwever.  we  have  found  the  Ca3- 
ore  agoatore  to  be  difTcrent  when  damage  is  cansad 
tingle  high  power  tnioowave  pulses  lasting  for  a  time  short 
compared  to  the  ihexmal  transfer  time  of  a  MESFET  gate^. 
In  1986  Anderson^  described  a  single  pulse  feilme  mecha¬ 
nism  called '‘snbsorfece  barnooi.’'  Be  based  his  arguments 
on  postmen em  observattons  at  the  vnr,H<y<  of  flje  damageo 
devices. 

Recnfi-3lioQ  feilnie  differs  sigmfk:ant1y  firom 
snhsurface  btrnooL  We  have  ifhinrtfiftrf  recuficaion  farime 
with  a  Qfnnhmflrtnn  of  diagnostir  techniques  never  before 
bcougfu  to  wear:  High*^>cedpitttogiapaj  was  ossd  to  clock 
tte  emigaon  of  cpdcal  tadiapoa  after  die  arrival  of  the 
bi^  power  microwave  poise,  cptical  spectroscopy  was 
employed  to  determine  the  source  of  the  optical  emissirms, 
and  rf  ^lectiosDopy  was  applied  to  monitor  the  genetadoa 
of  banuooics  of  the  incident  radiation  as  a  fimetion  of  time 
dmiug  each  test  pulse.  Devices  tested  inclnded  packaged 
low-noise  gtnpUfiers,  and  power  amt^ifiers. 

The  high-^ieed  pbotogra^ms  were  obtmoed  with 
an  optical  r^eak  camera.  Flgme  4  b  a  supexpo^on  of  the 
image  of  me  ciicalt  with  (he  streaked  Image  of  a  ^laik  that 
bomed  at  the  point  martrai  with  the  mrerw.  The  fignte 
shows  that  the  spark  first  became  vi^le  aLotn  ps 
alter  the  amval  of  the  miciowave  poke.  This  is  somewhat 
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Soak  photograph  showing  the  time  development  of  a 
qierk  00  a  MESFET.  The  spark  lasted  looger  than  the 
if  tea  pulse  that  aipposedly  gtacraied  the  qaik. 


tsmarfcable  coasideang  that  the  rniaowa^'e  poise  only  lasted 
for  80  ns.  Even  more  revealing  is  the  growth  is  hrigfame^ 
of  the  qffidc  for  the  remahnng  6  ps  trf  the  sneak  camfa-a 
sweep  period.  TbeciicQiishowrinFgure4was?power 
amplifier  fnEy  biased  for  nocnal  operation  KectificatioQ 
fedure  occcced  possihly  through  if  heating  of  the  Schooky 
juDCticsL  Afier  die  if  pulse  had  passed,  the  itvsise  bias  po- 
lenhab  were  able  to  drive  ctments  diat  gradually  heated 
the  gats  GO  mcandcxcocs.  This  mcandesceoce  was  re¬ 
coded  by  d"^  soeak  camera  as  a  ^adc  of  gtowsig  intsu- 
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Nitroges  lioe  nwfiarinn  oniiied  by  a  qwrk  fiom  an  un¬ 
biased  device  mdicturs  an  arc  between  elecaodes  on  the 
surface  of  the  MMIC. 


A  spearograph  of  the  Eghi  emiaed  fiom  a  similar 
device,  hisewise  biased  for  nocmal  t^iaaQOD.  mdicaied  a 
ccciiunous  9ecm>m  over  the  range  of  the  field  of  view 
(432  <  wavelength  (X)  <  632  nanometers^  Rgureodiows 

spectrom  obtained  fiom  an  unbiased  cL  Ibe  peaks  in 
the  scan  of  mnaisny  vs  waveleogib  are  etwyd  by  band 
emission  foam  fatciicd  nibogea  molecoles.  The  positioas 
of  (he  known  bands  are  idendfied  in  the  fignte.  The  dose 
agteeoieoi  between  the  observed  peaks  and  the  known 
jiCsinonsoftberunpgeD  bands  is  evidence  dot  atring  occurs 
between  electrodes  on  the  tBpface  of  die  chip  A  quod- 
graph  of  the  damaged  legiao  showed  signs  of  mild  heating 
but  no  anode-cathode  pits  ootmrjoly  «  «rv-»at>vt  with  arcs 
between  decttodcs  in  air.  We  somused  that  if  jonle  heal¬ 
ing  of  the  mjBfa-tfil  in  the  hi^i-fidd  lepon  of  the  gate  neck 
eocooraged  field  emf«a3An  of  elecboos  lasting  for  some 
fracrinn  of  the  if  poise. 

Hgure  6  diows  the  history  of  the  teedfied  second 
harmonic  emined  firmn  the  it^u ;  of  the  device  on  two 
SDCcesave  shots:  the  poise  neoeu.iig  the  damaging  pulse. 
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Secood  harmnnk  producaoc  during  the  damagning  test 
and  during  the  prergditig  test.  The  vanishing  of  toe  second 
ttyrmonir  is  IW  SiglEUDie  of  n-rrifirarinn  faflarc. 

and  the  damaging  poise  itself.  Ptior  to  damage,  the  ampU* 
tnde  of  the  bannonks  are  propartioaal  to  that  of  the 
damental.  but  cm  the  damaging  shot,  the  bannomes  vanish 
ead;  wbOe  die  ftmdamenal  cqntinoes  cm.  This  is  the 
ggnamne  of  recoficaium  fftilntr..  Hanoonics  are  geneiaied 
by  the  dippang  action  of  the  gaiejuncdoo  which  acts  as  a 
diode  alternate y  dtivea  into  fbrwrard  conduction  doting  die 
torward  half  cyck  and  into  avalanche  breakdown  during 
die  levetse  half  cycle.  When  the  JancDon  ceases  to  aa  as 
a  diode,  the  dipping  aedao  ceases  and  the  hannooics  vani^ 
Miaosccpically.  tte  cause  of  lectificaiicia  &iluie  is  saH 
unknown.  The  jonednn  metalliirsy  may  be  aUeted  by 
bearing  or  the  pheaomcaoa  may  be  leiaied  to  a  type  of 
CDtreot  mode  second  breakdown  observed  by  Van  Lint  et 
af.  This  year  we  have  plans  to  study  both  possibie  causes 
the  metallurgy,  nsing  tunnel iug  det^n  spectroscopy,  and 
the  breakdown,  using  broad  band  video  pulses  in  place  of 
die  if  poises  used  in  the  pR^vioos  ezpentnetus. 

While  tbese  smdies  arc  in  progres,  we  ate  ei- 


ploimg  the  knowiedgewe  have  gained  on  damage  mecha¬ 
nisms  by  having  some  off-the-shelf  MMIC  amplifiers  al¬ 
tered  to  withstand  greater  mcideiit  energy.  This  first  gen- 
eiatioo  of  baideoed  devices  is  srbcduled  for  testiog  in  July 
C^1991. 
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Rectification  failure  has  been  found  to  be  a  major  high-power  microwave  fail¬ 
ure  mechanism  in  amplifiers  employing  gallium  arsenide  metal  field  effect  transistors 
(GaAs  MESFETs).  The  mechanism  results  in  the  destruction  of  the  rectifying  prop¬ 
erty  of  the  Schottky  junctions.  A  junction  that  has  so  failed  can  no  longer  block  the 
flow  of  current  in  either  direction  and  acts  as  an  ordinary  ohmic  resistor.  Catastrophic 
failure  follows  from  current  driven  by  bias  voltages  through  the  failed  junction. 

Previous  investigators  employing  trains  of  lower  power  pulses  had  identified 
metal  migration  as  a  key  failure  mode  for  extended  exposure  to  moderate  levels  of 
microwave  radiation  (1-5).  However,  we  have  found  the  failure  signature  to  be  different 
when  damage  is  caused  by  single  high  power  microwave  pulses  lasting  for  a  time  short 
compared  to  the  thermal  transfer  time  of  a  MESFET  gate  (6).  In  1986  Anderson  (7) 
described  a  single  pulse  failure  mechanism  called  “subsurface  burnout.”  He  based  his 
arguments  on  postmortem  observations  of  the  surfaces  of  the  damaged  devices. 

Rectification  failure  differs  significantly  from  subsurface  burnout.  We  have 
identified  rectification  failure  with  a  combination  of  diagnostic  techniques  never  before 
brought  to  bear.  High-speed  photography  was  used  to  clock  the  emission  of  optical 
radiation  after  the  arrival  of  the  high  power  microwave  pulse,  optical  spectroscopy  was 
employed  to  determine  the  source  of  the  optical  emissions,  and  rf  spectroscopy  was 
applied  to  monitor  the  generation  of  hcirmonics  of  the  incident  radiation  as  a  function 
of  time  during  each  test  pulse.  Devices  tested  included  packaged  low-noise  amplifiers, 
«,nd  power  amplifiers.  The  tests  proceeded  by  the  direct  injection  of  a  single  pulse 
of  microwave  radiation  into  the  input  port  of  the  device.  The  pulses  typically  lasted 
for  100  ns.  Testing  would  start  at  power  levels  too  low  to  inflict  damage  and  would 
proceed  with  successively  higher  power  pulses  until  damage  was  observed.  Devices 
were  tested  both  with  and  without  bias  voltages  applied. 
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ABSTRACT 

Simulations  of  n-channel  silicon  MOSFET’s  with  gate 
□gths  between  0.05  and  0.25  microns,  performed  using  a  two 
mensiona]  self-consistent  energy  transport  model,  show  that 
ith  constant-field  scaling  maximum  electron  temperature  de- 
eases  significantly  as  channel  length  decreases  for  channel 
ngths  less  than  0.18  microns.  This  result,  indicates  that 
lort-channel  devices  can  be  more  reliable  than  similar  long 
>annel  devices. 

INTRODUCTION 

Reduction  in  MOSFET  channel  lengths  and  consequent 
icreases  in  peak  electric  fields  have  led  to  reductions  in  de¬ 
ice  reliability  through  hot-electron  effects  such  as  threshold 
oltage  shift  and  parasitic  gate  current.  However,  it  is  im- 
lortant  to  realire  increases  in  electron  energy  require  not  only 
hat  high  electric  fields  be  present,  but  also  that  sufficient  scat- 
ering  events  occur  to  randomize  electron  momentum  and/or 
nergy.  As  channel  lengths  are  shortened,  chances  that  such 
nteractions  are  “sufficient”  arc  reduced,  and  for  a  given  field, 
tlectrons  in  a  shorter  de\ice  may  acquire  less  energj’-i.e.  be- 
x)me  less  “hot”-  than  electrons  in  a  longer  device.  In  fact,  re- 
*nt  experimental  results  indicate'  that  hot-electron  induced 
lubstrate  current  in  Si  MOSFET’s  is  reduced  with  channel 
lengths  less  thi.n  0.15  micron.  Therefore,  the  hot-electron  reli¬ 
ability  problem  may  not  monotonically  become  more  severe  as 
channel  lengths  are  decreased. 

In  order  to  verify  this  apparent  shorter-channel  benefit, 
snd  to  quantify  where  beneficial  results  might  be  expected,  we 
have  investigated  the  behavior  of  n-channel  silicon  MOSFET's 
■rith  gate  lengths  under  0.25  micron,  using  techniques  which 
yield  true  electron  energies.  These  techniques  are  very  efficient, 
consuming  only  a  few  minutes  of  SUN  3/60  epu  time  per  device 
bias  point.  Tbiis,  submicron  device  reliability  predictions  can 
be  made  using  an  engineering  workstation. 

THEORETICAL  TECHNIQUES 

N-channel  siUcon  MOSFET’s  with  channel  lengths 
between  0.05  and  0.25  micron  were  simulated  using  the  pro¬ 
gram  UMDFET.*’*  UMDFET  self-consistently  solves  the  Pois¬ 
son,  continuity,  and  energy  bsdance^  equations  in  two  dimen¬ 
sions.  Changes  in  electron  energy  due  to  acoustic  and  optical 
phemoa  scattering,  heat  flow,  and  diffusion  are  included. 

The  devices  in  our  computer  experiment  were  designed  so 
■s  to  keep  the  magnitude  and  shape  of  the  lateral  electric  field 
the  same  in  all.  Doping  concentrations  and  bias  voltages  on 
the  gate  and  drain  were  therefore  adjusted,  as  shown  in  Table 
1.  A  mavimiim  lateral  electric  field  of  200K V / cm  was  imposed 
on  all  devices.  The  gate  and  drain  bias  voltages  for  each  device 
were  kept  almost  equal  so  that  the  transverse  electric  field  was 
negligible  in  comparison  to  the  lateral  electric  field  near  the 
drain,  where  the  electrons  reach  their  maximum  temperature. 
The  ambient  temperature  for  all  of  the  devices  simulated  was 
300K. 


RESULTS 

Sample  profiles  for  lateral  electric  field  at  the  surfaces  of 
three  of  the  simulated  devices,  and  the  resulting  curves  of  sur¬ 
face  electron  temperature  as  a  functon  of  distance  from  the 
drawn  gate  edge,  are  shown  in  Figure  1.  The  results  are 
summarized  in  Figure  2.  These  results  show  that  for  channel 
lengths  below  0.18  micron  the  maximum  electron  temperature, 
related  to  the  maximum  energy  a'vailable  to  cause  reliability 
problems,  decreases  as  channel  length  is  decreased.  For  chan¬ 
nel  lengths  above  0.18  micron  the  maximu:  ■  ,-ctron  energy 
changes  little.  For  very  long  channels  elec..jn  temperature 
is  expected  to  approach  a  maximum  value  around  4300K,  the 
temperature  reached  in  a  constant  200KV/cm  field.® 

DISCUSSION 

The  electron  temperature  shown  in  Figures  1  and  2  is  an 
average  quantity  which  is  a  measure  of  the  “randomness”  of  the 
electron  energy,  i.e.,  the  spread  of  electron  energies  around  an 
average  value.  Thermalization,  an  increase  in  this  randomness, 
is  caused  by  scattering  events  as  the  electrons  travel  through 
the  device.  A  higher  electron  temperature  at  any  point  implies 
that  more  electrons  have  undergone  randomizing  collisions  up 
to  that  point.  Some  of  the  electrons  may  have,  on  balance,  ac¬ 
quired  energy  from  the  lattice;  some  may  have  lost  it.  Electrons 
with  energies  larger  than  the  lattice  temperature  are  “hot”.  If 
hot  enough,  these  electrons  cause  reliability  problems.  In  any 
case,  the  average  electron  temperature  (plotted  in  Figs.  1  and 
2)  is  an  indicator  of  the  spread  of  energies  to  be  found  in  the 
electrons;  the  larger  this  spread,  the  greater  the  number  of  hot 
electrons. 

In  a  MOSFET  electrons  are  accelerated  by  the  lateral  elec¬ 
tric  field  over  the  length  of  the  channel.  If  the  channel  is  very 
short,  the  field  has  only  a  limited  time  to  act  on  the  electrons, 
and  the  electrons  can  suffer  only  a  relatively  few  randomizing 
collisions.  However,  as  channel  length  is  increased,  the  elec¬ 
trons  have  time  to  reach  a  maximum  value  of  energy  (i.e.,  the 
range  of  energies  which  they  exhibit  reaches  some  equilibritim 
value).  Figure  2  shows  that  this  critical  length  is  around  0.18 
micron  for  the  field  profiles  we  have  used.  This  is  also  the 
value  of  channel  length  at  which  velocity  overshoot,  a  related 
effect,  has  been  shown  to  appear  in  short-channel  MOSFET’s.* 
A  more  complete  discussion  of  the  relationship  of  electron  tem¬ 
perature  to  device  reliability  has  been  given  elsewhere^. 

CONCLUSIONS 

In  very  short-channel  devices,  the  transit  time  of  electrons 
through  the  high  electric  field  may  be  shorter  than  the  en- 
ergy  relaxation  time.  When  this  happens,  the  electrons  do 
not  have  enough  time  to  fully  exchange  energy  with  the  lat¬ 
tice.  The  result  is  a  lower  average  electron  temperature  than 
would  occur  in  a  similar  long-channel  device.  Hot-electron  re¬ 
lated  MOSFET  reliability  problems  should  therefore  become 
less  threatening  for  silicon  MOSFET’s  with  very  short  chan¬ 
nel  lengths,  even  if  the  electric  fields  and  drain-source  voltages 
in  these  devices  remain  relatively  high.  For  devices  scaled  as 
those  here,  improvements  should  be  noticed  for  chaimel  lengths 
under  about  0.18  micron. 


We  would  like  to  thank  T.  Urai  for  his  help  with  UMDFET 
and  many  discussions  on  electron  transport  phenomena. 
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Figure  1.  Lateral  electric  field  (long  dash  short  dash)  and  elec¬ 
tron  temperature  (solid)  profiles  vs.  distance  from  drawn  gate 
edge  for  (a)  0.05  micron  gate  length,  (b)  0.15  micron  gate 
length,  and  (c)  0.25  micron  gate  length  n-channel  silicon  MOS¬ 
FET’s. 
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Nd=max.  doping  concentration  at  drain  and  source  (cm-3) 
Na=implant  concentration  in  channel  (cm-3) 
xj=drain  and  soiuce  junction  depths  (micron) 
tox=oxide  thickness  (nm) 

Vgsgate  bias  (V) 

Vdsdrain  bias  (V) 

Vf=Vb=0.0,  where  Vs  and  Vb  are  the  bias  voltages  of  the 
■ource  and  bulk  (Device  doping  profiles  obtained  using 
process  simulator  from  MINIMOS.  Drain  and  source 
junctions  are  diffused  As,  chazmels  are  B  ion  implants) 


Table  1.  Device  parameters  of  the  simulated  MOSFET’s. 
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chsnntl  length  (alcron) 

Figure  2.  Maximum  electron  temperature  at  surface  of  chan¬ 
nel  vs.  channel  length  for  0.05,  0.08,  0.1,  0.15,  0.18,  and  0.25 
gate  length  n-chaimel  silicon  MOSFET’s  with  peak  lateral  elec¬ 
tric  field  of  200KV/cm.  Error  bars  of  -f/-150K  reflect  under¬ 
certainties  in  scattering  parameters  and  simulation  program 
convergence  tolerances. 
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A  method  for  calculating  impact  ionization  coefTicienls  by  solving  the  Boltzmann  transport 
equation  is  presented.  The  distribution  function  is  taken  to  be  expressible  as  a  Legendre 
polynomial  expansion,  which  is  substituted  into  a  Boltzmann  equation  that  incorporates  the 
effects  of  nonparabolic  band  structure,  deformation-potential  phonon  scattering,  and  impact 
ionization.  The  resulting  Boltzmann  equation  can  be  expressed  m  a  linear  form,  and  solved 
using  sparse-matrix  difference-differential  methods.  Ionization  coefficients  are  obtained 
directly  from  the  distribution  function.  Calculated  values  for  the  ionization  coefficients  agree 
very  well  with  experiment  for  electrons  in  silicon 
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RELY,  a  new  MOSFET  tunulaior  for  predicting  hoi 
electron  induced  degradatioo  problems  is  presented.  RELY 
uses  the  distribubon  funebon  emulated  by  solving  the  Energy 
and  Boltzmann  transpon  equabons  id  allow  predicbon  of  the 
bme  dependence  gate-leakage  current,  oxide  charge 
deposition,  and  device  I-V  charaoerisbes.  Results  agree  with 
experiment  and  Monte  Carlo  calculations,  and  are  o^ned  in 
only  several  CPU  minutes. 


INTEQDLCT1Q.N 

We  present  here  RELY,  a  novel,  physics-based  CAD  tool 
for  predicting  MOSFET  leliabiliiy  ptoblerrts  As  device 
dimensions  conbnue  to  shrink,  many  forms  of  device 
degradabon  will  be  prevalenL  Cate-leakage  current,  threshold 
voluge  shifts,  and  transcortductance  thifu,  arising  from  hot- 
electron  damage,  are  parbcularly  threatening,  these  phenomena 
degrade  device  operating  characterisbes  over  broe,  and  may 
ultimately  kad  to  device  failure.  Models  previously  proposed 
for  predicting  teliabibry  problems,  including  the  popular 
Richardson's  Equabon  Method  for  emission  over  a  barria,  are 
often  based  on  the  incorrea  assumpbon  that  die  electron 
^stribubon  funebon  is  Maxwellian  (]}.  Other  methods, 
including  Monte  Carlo  simulabon,  are  more  physically 
rigorous,  but  require  prohibibvc  amounts  of  CPU  bme,  and  are 
thus  inappropriate  for  use  in  most  CAD  tools  for  prediebng 
device  reliability  |2] 

The  new  device  reliability  simulaicr,  RELY,  is  both 
physics- based,  and  oompuubonally  effrcieni.  RELY  uses  a 
new.  Energy  Transpon-Boltzmann  Transpon  Equation  (ET- 
BTE)  model  Since  RH.Y  is  physics-based,  the  ne^  for  many 
fitbng  parameters,  often  associated  with  predicting  device 
degr^bon,  is  eliminaied.  ‘The  ET-BTE  model  uses  device 
physics  -  icanering  mecharusms  and  oonparaboUc  band 
structure  -  to  calculue  the  average  electron  energy  and  the 
electron  distribubon  function  throughout  a  device.  By  utilizing 
this  informabon,  RELY  predicts  MOSFET  ^iie-leakage  current 
and  threshold  voluge  shifi  lo  ascertain  dieir  effects  on  device 
operabon  as  a  function  of  atressina  time  Excellent 
agreeineni  with  experiment  is  anaintd;  calcuiabons  have  also 
bttn  shown  lo  agree  with  Monte  Carlo  stmulabons  while 
requiring  less  than  T/lOO  the  CPU  time  to  evaluate. 

MODEL  DESCRIPTION 


ET.BTE  Model 


The  ET-BTE  model  consists  of  five  equabons:  the  Energy 
Transpon  equabon  (ET)  (1),  and  the  Boltzmann  Transpon 
equabon  (B'lE)  (2).  as  well  as  the  Poisson  (3),  conbnuity  (4) 
and  current  (S)  equabons. 
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where  yr  is  the  elecirostabc  potential,  n  and  p  are  electron  and 
hole  concentrations,  N  is  the  net  fixed  charge  density,  R  is  the 
net  recombinabon  rate,  is  the  electron  current  density  >  - 
enu).  V.  aa,  E,  f(k),  and  k,  are  the  electron  average 

velocity,  electron  average  energy,  electric  field,  distribubon 
funebon,  effecbve  electric  field,  energy  lelaxabon  time  and 
wave  vector  tespeebvely.  The  subscripts  or  and  iti  in  the  BTE 
icprcseni  the  effects  of  acoustic  and  intcrvalley  phonon 
acanering  on  the  distribubon  funebon 

To  fivnliute  predicting  nme-dependem  device  reliability 
problems,  RELY  calculates  the  change  in  currents  and 
elecirostabc  potential  due  to  boi-eiectron  degradabon  as  a 
function  of  operation  nme.  To  begin  these  calculations,  the 
inibal  values  (bme  M))  for  dectrosabc  poienbaJ  and  earner 
concentrabons  arc  provided  a  2-D  Poisson-continuity 
equabon  solver,  such  as  MINIMC^  [3].  With  concentrabon  and 
poeenbal  as  input,  RELY  then  aolves  die  energy  balance 
equabon  (1)  in  two  dimensions  To  solve  equabon  (1),  a  finite 
difference  (FD)  scheme  is  ^iplied  Discrebzabon  is  performed 
at  boundary  elemenu  by  modification  of  the  FD  formulae  using 
die  muTor  image  lechn^ue  [4].  This  lechi^ue  provides 
accurate  use  of  equation  (1)  ^  elimin^g  the  typical 
lequiiement  that  carrier  and  lanioe  tcmpeiatures  be  equal  at 
boundaries.  Finally,  a  constant  value  for  energy  relaxatian 

time,  shown  by  Monte  Garb  calcuiabons  to  be  2.6xl(r^^  see, 
is  used  |5].  The  aolubon  of  equabon  (1)  provides  the  ewage 
electron  energy  as  a  funebon  of  the  device's  space  coordiMtes 
Device  regions  of  high  electron  energy  arc  potential  problem 
veas. 

To  find  the  duoibubon  fkmebon  at  ea^  channel 
coordinate,  die  average  energy  found  above  is  naed  in 
conjunction  with  aotuboru  lo  die  homogeneous-field  Boltznann 
Tran^nn  Equabon  (2).  Due  id  its  coorplexity,  we  present  here 
only  a  symbolic  icprejuiation  of  dr  BTE.  Within  RELY,  die 
hotnogeneous-field  BTE  ia  forroulated  using  Legendre 
poiynotnials,  while  for  the  effects  of  silicon’s  i 


nonpanboiic  band  nructure,  iniervalley  and  acoustic  phonon 
icanavig  (6]  The  details  of  the  phonon  collision  intepils  arc 
expressed  in  lenns  of  deformation  p«entii]  scattering  theory 
V^ues  for  deformation  potentials  arc  identical  to  those  used  for 
Monte  Carlo  calculations  [7]  To  account  for  space  depeni^ce 
an  effective  electric  field  is  used  as  input  to  equation  (2). 

is  obtained  by  mapping  position-dependent  avenge  energy, 
found  from  the  solution  to  equation  (1),  to  the  correqxnding 
eiectnc  field  on  die  homogeneous-average -energy  versus 
clecov-field  curve.  This  curve  was  calculated  by  Monte  Carlo 
analysis  beforehand.  After  inputting  E^  RELY  numerically 

solves  the  BTE  using  the  finite  difference  method 

The  solutions  to  Ae  above  five  equabons  provide  the 
avenge  energy  and  the  electron  disiribubon  function  throughout 
the  device,  as  well  as  the  I-V  characterisbes 


MOSFET  gale  leakage  current  appears  when  electrons 
gain  enough  energy  to  surmount  or  tunnel  through  the  oxide 
barrier  The  magnitude  of  thermal] y-emined  gate-leakage  current 
is  tpven  by  the  number  of  electrons  which  are  energeuc  enough 
to  surmount  the  oxide  barria,  mulbplied  by  their  respecbve 
velocities  in  the  direebon  of  the  oxide  Mathemabcally  this  is 
given  as 

y,(x)*J  en(x)  o/k)F(x,k)d\  (6) 

where  y^fx)  is  the  gate  current  density,  n(x)  is  the  electron 
concentrabon.  k^  is  the  irbnimum  electron  wave  vector 

necessary  for  surmounting  the  oxide  barrier,  u^fk)  is  the 
electron  velocity  perpendicular  so  the  interface,  F(xjt)  is  the 
space -dependent  tnomenturti  distribubon  funebon. 

By  performing  a  great  deal  of  algebraic  manipuiabons. 
which  are  made  di^cult  by  silicon's  ellipsoidal,  nonparabolic 
band  structure,  we  conven  the  above  equation  so  an  explicit 
expression  for  calculabng  gate  current  density  in  MOSFETs 
[8]  This  new  expression  is  similar  to  Richardson's  equabon  for 
thermal  emission  over  a  barrier.  However,  unlike  Richardson's 
equabon  which  is  based  on  a  Maxwellian,  the  new  expression  is 
d^ved  from  the  physics-based  hot-electron  distribubon 
funebon,  which  is  found  from  the  ET-BTE  model  described 
above  This  new  expression  is 


of  chatvtel  ooonlinaie.  Total  gate-leakage  current  is  obtained 
by  integmbng  J^(t)  along  the  channel: 


f*«i 


‘  3h’ 


where  ti\  and  m|  are  silicon’s  transverse  and  longitudinal 
cfrecbve  masses  lespecbvcly.  7(c)  is  Eicon's  nonparaboUc 
dispersion  relabon  (7],  is  the  oxide  barrier  height.  f(x.c)  is 
the  ^ace-dependent  energy  distribubon  function,  which  has 
units  of  cm3;  ji  i},^  probability  of  finding  an  electron  in  the 

vicinity  of  energy  e  and  location  a. 

The  effects  of  Schooky  barrier  lowering  and  electron 
mnncling  are  accounted  for  by  adjusting  the  oxi^  barrier  height 
according  to  the  following  expteuion  [9j: 

♦p-3.I-b£^-aE^  (S) 

Eqj  is  the  electric  field  in  the  oxide.  1  and  b,  which  have  been 
dcterrauied  by  oouyiarison  arith  experiinent.  have  the  values 
1 X  KT*  e(CTr*v)'^  and  2.6x KT*  e(cmV)*^  respectively. 

R^Y  oumerically  iuegraies  equation  (7),  while 
inserting  proper  values  for  Rm)  determined  by  the  ET-BTE 
nKthod,  ID  asoenain  gate  leakage  civreni  density  as  a  function 


W|  is  the  width  of  the  MOSFET  pte.  x,  and  x^  are  the  source 
and  drain  coordinates  respectively. 


A  small  fraction  of  the  hot  electrons  which  surmount  the 
oxide  barrier  become  trapped  in  oxide  sates.  The  oxide  charge 
deposition  gives  rise  to  higha  channel  resistivity  and  increased 
threshold  voluge.  and  ultimately  affects  the  basic  I-V 
characteristics  of  the  device  RELY  predicts  time-dependent 
oxide  charge  buildup  by  using  the  foUowuig  equabon  which 
describes  the  exponential  lelabonship  between  the  occupation  of 
oxide  states,  the  gate  current  density,  and  the  device  soessing 
time  19]: 

i 

where  is  the  density  of  type  1  oxide  traps,  Af^(t)  is  '.he 
number  of  oxide  sates  filled  at  time  t,  is  the  cross  section  of 
oxide  trap  i 

Using  the  previously  calculated  gate  current  density, 
RELY  evaluates  equation  (10)  to  predia  values  for  oxide  charge 
accumulation  with  time.  The  default  values  for  trap  densities 
and  cross  sections  are  [10] 

O)  -  1.71x10-'’,  O2  •  8.81x10-'’,  03  *  1.93x10-'* 
^p,,«9.65xl0",  ■  9.43x10",  -  1.14x10” 

The  trapped  charp  distributions  in  the  oxide  ofen  cause  a 
localized  perturbtmon  on  the  potential  and  carrier  distributions 
near  the  d^  region  of  the  device.  This  perturbation  potential  is 
determined  in  two  dimensions  using  the  time  perturbation 

method  [11].  At  tiiDe  t|,  the  Poisson's  equation  can  be  wrinen 

as 

V  Vb  (x,y,i,)»p(x,y,t,)/c„,  ,  in  oxide  (11) 

V  Vd  (x,y,t|)»e(n'-p'-N)/t,  ,  in  silicon  (12) 

where  f  (x,y,t|)  is  the  new  potential  because  of  the  penuibation 

caused  by  the  napped  chvges,  a'  nd  p'  ate  ibe  new  tkxiraii 
and  bole  concentrations,  p(x,y,t|)  («  e  Af0j(x,y,t|))  is  the 
trapped  oxide  charge  rjimUtwi  above. 

The  new  potential  can  be  deumyosed  into  ■  time 
independent  pan  and  a  time  dependent  part  At  time  t|, 

♦  (*.y.«|)  -  9<*.y)  *  B(»-y.«,)  (13) 

where  u(x,y,t|)  is  the  pemrbatioo  potential  at  time  tj. 

We  nan  with  a  potential  disiribution  and  electron 
ooncentrations  obtained  from  MINIMOS.  Then,  the 
perturbation  potential  is  obtained  after  nibstituting  relation  (13) 
into  (11)  and  (12),  and  lineariang  the  resulting  equations  Akfro- 
die  perturbation  potential  is  determined,  local  correctiofu  are 
made  k>  ihe  carrier  conctnirations  In  linae,  at  mtae  charge 
becomes  sapped  in  die  eude,  we  penodically  use  ihc 
penurbation  method  m  ic-«valuate  dn  Axsson  and  continuity 
aquations.  At  each  iiention,  die  new  I-V  characteristics  and 
dueshold  voltage  shift  ac  ujng  updated  potential 

distributions  and  charge  oonoeaiiBtions.  The  Sow  diagTam  of 
RELY  ii  shown  in  Figure  1. 


Figure  1  The  Flo*  Ditgrun  of  RELY 


RESULTS  AND  DISCUSSIONS 

To  demonstrate  the  capabilities  of  RELY,  we  smuiaied 
n-channel  MOSFETs  on  a  SUN  desktop  computer  Fipre  2 
shows  energy  versus  position  along  the  channel,  areas  of  large 
electron  energy  indicate  repons  subject  lo  increased 
degradation.  Figure  3  compares  ^stribution  functions  calnilated 
by  RELY  (Legendre  poijmofnial  lech^ue),  the  Monte  Carlo 
method,  and  the  Maxwellian  distiibubon.  l^le  there  is  food 
agreement  between  Legendre  polyrxrnu^  imd  Monte  Carlo 
techniques,  the  Maxwellian  disagrees  significantly  ower  all 
enerpes  Such  agreement  between  Monte  Carlo  and  Legendre 
polynomial  techniques  indicates  that  the  new  ET-BTE  method 
does  indeed  predict  realistic  values  for  distribution  function. 
RELY  predicts  values  for  gate  current  which  match  very  well 
with  experiment  for  over  four  order  of  magnitude  as  shown  in 
Figure  4.  Unlike  other  models,  which  ignore  the  actual  physi^ 
mechanisms.  RELY  attains  agreement  for  a  wide  tinge  of  bias 
conditions,  using  the  poorly-known  phonon  defonnation 
potential  as  hs  only  adjustable  parameter.  When  more  reliable 
values  for  these  phonon -electron  coupling  factors  are  available, 
we  are  likely  to  raider  unnecessary  the  one  fining  parameier 
currently  required  by  RELY  for  gate  current  predictions,  ngure 
S  compares  oxide  charging  of  the  test  device  after  eperabon 

times  of  0.01,  1.0,  and  100  hours.  Reasonable  values  for 
Kinporal  evolution  of  drain  current  luve  also  been  obtained  as 
shown  in  Figure  6.  The  presented  method  of  oxide  charging  has 
also  been  ^iplied  to  timulate  EPROM  programming 
chsractoistics  a^  exccUeni  agreement  with  caperiment  has 
been  obtained  u  shown  in  Figure  7(12]. 


Figure  2  Average  energy  versus  position  along  the  dhannel  The 
device  is  biased  at  V^-6.0V  and  Vjj*7i)V.  LAV  > 

L2titD  /  S0.8)jm,  tgi  w  150A.  Gate  oxitfc  begins  at 


Figure  3  The  Le^ndre  jKilynomia]  distribuiion  fimetion  and 
Maxwellian  di^bubon  ve  compared  with  the 
distribution  funebon  calculated  using  liooie  Carlo 
•chnique.  All  distributions  are  for  an  avoage  energy 
of0.72ev. 
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houn  of  operation, 
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Figure  6  Calculated  drain  current  demdaiion  (dl  /  I)  versus 
tocssing  time.  Stress  condition  Is  Vds»8,0V  ar.d 
V|s-3.0V.  "nie  characteriiation  is  with  Vds»0.1V  and 
. Vgs«3.0V.  L/W  ■O.^im/SO.Sum,  tp,  ■  ISOA. 
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CONCI 

I^Y.  a  Cad  tool  which  combines  phy 
CTgtneer^  pmcncaliiy.  has  been  develotSd 
lanuies  of  calculation  on  a  wort-station  cot^te 
^'•fcpenda.t.  ta^leciron  wSed 
W^len^  as  well  as  mean  time  so  failure,  in  de 
fabnea^  SiKh  nmulations  riiould  prove  va 
deaign  of  reliable  lubmicrao  aemicooducuir  device 
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abstract 

This  paper  describes  a  novel  technique,  which 
combines  the  attributes  of  the  Energy  Transport  and 
Monte  Carlo  methods,  for  determining  MOSFET  gate 
current  that  arises  from  electron  heating  in  the  device. 
This  method  is  based  upon  a  non-Maxwellian  hot- 
electron  distribution  function  found  from  Monte  Carlo 
simulation,  and  utilizes  a  physically-calculated  average 
electron  energy.  Calculated  value's  for  gate  current  of  a 
sample  subnticron  MOSFET  device  show  reasonably 
good  agreement  with  experiment  without  the  need  for 
any  fitnng  parameters. 


INTRODUCTION 

Hot  electron  injection  into  the  gate  oxide  is  the 
underlying  mechanism  largely  responsible  for  MOSFET 
gale  leakage  current,  MOSEFTT  performance  degrada¬ 
tion,  and  EPROM  programming  To  accurately  investi¬ 
gate  these  hot-electron  phenomena,  the  energy  distribu¬ 
tion  function  is  nevded.  However,  neither  the  conven- 
bonal  Drift-Diffusion  (DD)  model  nor  the  more  complex 
Hydrodynamic  Energy  Transport  (ET)  model  provide 
the  distribution  function.  Therefore,  much  effort  has 
been  invaied  in  predicting  oxide  injection.  Many 
models  previously  proposed,  including  tJie  popular 
Richardson's  Equation  method  [1,2],  are  based  on  the 
disputed  assumption  that  the  distribution  function  is 
Maxwellian.  Other  approaches  do  nut  make  such  an 
assumption,  but  are  computationally  too  expensive  to 
be  of  day-to-day  use  for  device  design  (3| 

In  this  paper  we  present  a  novel  technique, 
which  couples  the  Energy  Transport  and  the  Monte 
Carlo  (MC)  methods,  to  determine  oxide  injection  and 
MOSFET  gate  current.  The  merits  of  this  new  Energy 
Transport-Monte  Carlo  (ET-MC)  approach  are;  (1)  it 
incorporates  the  semiconductor  band  structure  and 
scattering  mechanisms;  (2)  it  renjuires  less  Cl'U  time 
than  a  standard  MC  analysis;  (3)  it  uses  a  robust  distri¬ 
bution  function  tail  population-enhancement  scheme; 
(4)  it  is  a  candidate  for  CAD  applications. 


*  ET-MC  MATHEMATICAL  MODEL 
The  ET-MC  is  composed  of  a  homogeneous-field 
Monte  Carlo  calculation  in  conjunction  with  the  Poisson 
(1).  continuity  (2).  current  (3).  and  energy  balance  (4) 
equations. 

+  —Ip  -  n  +  N)  =  0  (]) 

=  fR  (2) 
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where  n,  p,  v,  (u,  E,  R,  and  are  the  electron 
concentration,  hole  concentration,  electron  cturent  den¬ 
sity.  electron  drift  velocity,  electron  average  energy, 
elKtric  field,  recombirution  rate,  mobility  and  diffusion 
coefficient  respectively. 


ET-MC  ALGORITHM 

To  evaluate  the  mathematical  model  described 
above,  we  perform  the  following  operations. 

1.  The  electrostatic  potential  and  carrier  concentrations 
within  the  device  are  calculated  by  a  2-D  Poisson- 
continuity  equations  solver,  such  as  MINIMOS  or 
UMDFETJ4,5J. 

2.  The  average  electron  energy  as  a  function  of  device's 
space  coordinates  is  obtained  by  solving  energy  balance 
equation  (4)  using  the  algorithm  described  in  reference 
16j.  The  electric  field  obtained  from  the  above  Poisson- 
continuity  solver  is  used  as  input  to  equation  (4). 

3.  An  effeciive-elecirie  fUld  profile  is  then  determined  to 
facilitate  determining  the  space-dependent  distribution 
function. 

4.  The  space-dependent  disiribubon  function  is  ascer¬ 
tained  using  homogencous-rield  MC  simulations  with 
the  effecbve  Field  as  input. 

5.  The  gate  current  density  is  found  by  employing  ther¬ 
mionic  emission  arguments  while  using  the  newly- 
calculated  distribution  function. 

The  novelty  of  this  work  is  contained  in  steps  3 
through  5  which  are  described  in  detail  below. 


Effective-Electric  Field  Profile 

To  account  for  the  effects  of  nonlocal 
phenomena  we  define  an  effective-electric  field  profile. 
We  determine  the  effective-electric  field  by  mapping  the 
average  energy  to  the  corresponding  electric  field  on  the 
homogeneous-electric  field  versus  average  energy  curve, 
which  is  shown  in  Figure  1.  After  performing  this 
operauon  for  each  space  coordinate,  an  effective  field 
profile  is  obtained.  The  space-dependent  average  energy 
has  already  been  found  from  the  solution  of  energy  bal¬ 
ance  equation  (4),  while  the  homogeneous-electric  field 
versus  average  energy  curve  was  calculated  by  MC 
analysis  beforehand.  The  effective  and  the  actual  electric 
field  profiles  in  the  channel,  parallel  to  the  interface,  are 
shown  in  Figure  2. 

Space-Dependent  Distribution  Function 

for  phenomena  which  occur  in  silicon  on  the 
time  scale  oi  die  energy  relavaiion  tune,  the  average 
eleciron  energy  is  an  accurate  yvarameter  to  use  as  a 
guideline  for  interpolating  homogeneous-field 
phenomena  to  space-dependent  phenomena  17).  We 
utilize  this  concept  to  deierrnine  the  space-dependent 
distribution  function  from  homogeneous-field  MC  cal¬ 
culations  More  specifically,  we  first  determine  points 
where  average  energy  is  greatest.  We  then  perform 
homogeneous-field  MC  simulations  at  each  of  those 
points  using  the  effective-electric  field  values  as  input. 

Calculating  the  High-Energy  Tail  of  the  Distribution 
Function 

To  obtain  the  tail  of  the  distribution  function, 
where  the  occupation  probability  is  very  small,  we 
divide  our  energy  domain  into  high  and  low  regions. 
We  first  generate  a  distribution  in  the  low  energy  region 
where  eleciron  states  are  easily  accessed.  Once  this 
region  is  statistically  well  defined,  i.e.,  enough  events 
have  occurred  to  make  the  curve  that  is  well  defined,  or 
smooth  around  this  point,  we  need  not  generate  any 

more  events  in  this  region,  it  is  considered  known. 
Now  we  start  the  electrons  at  the  next  energy  region. 
This  region  begins  at  the  highest  energy  point  where 
this  initial  distribution  function  is  well  defined.  We  let 
this  stage  of  the  Monte  Carlo  proceed  as  usual,  but 
when  the  electron  energy  falls  below  the  energy  break¬ 
point  we  have  picked,  it  is  thrown  away  and  another 
electron  is  started  again  at  this  breakpoint.  The  pari  of 
the  distribution  function  generated  by  this  sequence  of 
events  must  now  be  properly  weighted  with  respect  to 
the  initial  distribution  function.  This  is  easily  done, 
because  we  already  know  the  'correef  probability  of  an 
occurrence  at  the  energy  breakpoint  from  the  initial  dis¬ 
tribution  function.  Thus,  we  match  values  /or  the  dis¬ 
tribution  function  at  the  breakpoint.  This  matching  is 
accomplished  by  multiplying  all  of  the  probabilities  in 
the  low  probability  region  by  a  single  weighting  factor. 
We  can  extend  the  distribution  function  to  higher  and 


higher  energies,  or  lower  and  lower  probability  regions, 
by  iterating  this  scheme  many  times. 

This  process  may  be  more  easily  understood  in 
mathematical  terms.  1/  Nj(c)  is  the  occupation  number 
in  the  initial  easily-generated  energy  histogram  at  an 
energy  e,  N-,(e)  is  the  number  obtained  at  energy  e  in 
the  tail,  and  e-  is  an  energy  breakpoint,  then  the  weight¬ 
ing  at  each  energy  breakpoint  ii  p»cr::xrzers  af  iolkiws 


1.  For  £  2: 

N(£)=Nj(£kN2tE)  15) 

2.  For  £  i  multiply  each  such  MU)  by  the  weighting 
factor 


N,(£,) 

Me  )  - - — - 

'  N,(£,}  M2U,) 


(6) 


3.  The  uiinornialized  distribution  function  is  then  given 
by 

fU)  =  N,U)  £:Si;  t7j) 


f(E)  =  N(e)W(£)  £  2  c,  (7b) 

Cate- Leakage  Current  Calculation 

MOSFET  gate  leakage  current  appears  when 
electrons  gain  enough  energy  to  surmount  or  tunnel 
through  the  oxide  barrier.  The  magnitude  of  thermally- 
emitted  gate-leakage  current  is  given  by  the  number  of 
electrons  which  are  energetic  enough  to  surmount  the 
oxide  barrier,  multiplied  by  their  respective  velocities  in 
the  direction  of  the  oxide.  Mathematically  this  is  given 
as 


;^U)  =  \enU)v^[k)ru.k)jh  (8) 

•« 

where  /^(z)  is  the  gate  current  density,  n(r)  is  the  elec¬ 
tron  concentration,  k^  is  the  minimum  electron  wave 
vector  necessary  (or  surmounting  the  oxide  barrier, 
v^{k)  is  the  electron  velocity  perpendicular  to  the  inter¬ 
face,  Flxjc)  is  the  space-dependent  momentum  distribu¬ 
tion  function. 


By  performing  a  great  deal  of  algebraic  manipu¬ 
lations,  which  are  made  difficult  by  silicon's  elLpsoidai, 
nonparabolic  band  structure,  we  convert  the  above 
equation  to  an  explicit  expression  for  calculating  gate 
current  density  in  MOSFCTs  |8].  This  new  expression 
is  similar  to  Richardson's  equation  for  thermionic  emis¬ 
sion  over  a  barrier.  However,  unlike  Richardson's  equa¬ 
tion  which  is  based  on  a  Maxwellian,  the  new  expres¬ 
sion  is  derived  from  a  physics-based  hot-elecuon  distri¬ 
bution  /unction,  which  is  found  from  the  ET-MC  simu¬ 
lation  described  above.  This  new  expression  is 


/^Cx)= 


1 

m  1/2 

m, 

, 

where  m^  and  are  silicon's  transverse  and  longitudi¬ 
nal  effective  masses  respectively,  7(e)  is  silicon's  nonpar¬ 
abolic  dispersion  relation,  is  the  oxide  barrier  height. 
/<x,£)  is  the  space-dependent  energy  distribution  func¬ 
tion,  which  has  units  of  0/1^;  it  gives  the  probability  of 
finding  an  electron  in  the  vicinity  of  energy  £  and  loca¬ 
tion  X.  The  effects  of  Schottky  barrier  lowering  and  elec¬ 
tron  tunneling  are  accounted  for  by  adjusting  the  oxide 

barrier  height  according  to  an  empirically  determined 
algorithm  [9].  Values  for  l^(x)  are  determined  by 
integrating  ecjuaiion  19),  while  inserting  prop)er  values 
for  jit^)  determined  by  the  ET-MC  method.  Total 
gate- leakage  current  Ig  is  obtained  by  integrating  /^(j) 
along  the  channel. 

RESULTS  AND  DISCUSSION 

We  demonstrate  the  new  ET-MC  method  with  a 
calculation  of  gate  leakage  current  ansing  from  an  n- 
channel  MOSFET  with  0.9^m  effective  gate  length  using 
a  SUN  desktop  computer.  The  dispersion  relation  and 
scattering  coupling  factors  are  taken  from  reference 
(10].  figure  3  shows  energy  versus  position  along  the 
channel;  areas  of  large  electron  energy  indicate  regions 
subject  to  increased  degradation  The  distribution  func¬ 
tion  at  channel  coordinates  l.CWSpm  and  1.061(im, 
corresponding  to  effective-electric  field  values  of 
202k V/cm  and  237kV/cm  are  plotted  in  Figure  4.  The 
excellent  resolution  in  the  high-energy  region  of  the  dis¬ 
tribution  function,  demonstrated  by  the  smooth  curves 
in  Figure  4,  is  attributable  to  the  tail  population- 
enhancement  method  discussed  above.  The  electron 
energy  distribution  functions  are  muluplied  by  the 
density  of  states  and  normalized  to  unity.  Values  calcu¬ 
lated  for  gate  current  using  the  presented  ET-MC 
method  and  Richardson's  equation  are  compared  with 
experiment  as  shown  in  Figure  5.  Figure  5  shows  rea¬ 
sonably  good  agrv'cmcnt  with  experiment  is  obtained 
using  the  ET-MC  method  without  the  nxvd  for  any  fit¬ 
ting  parameters;  Richardson's  method  is  in  error  by 
approximately  orders  of  magnitude.  While  reason¬ 
ably  good  agreement  with  experiment  was  attained,  we 
believe  our  method  will  be  even  more  accurate  when  an 
even  more  detailed  MC  program  is  used.  Furthermore, 
calculation  times  can  be  significantly  reduced  by  pro¬ 
viding  a  database  of  distribution  functions,  each 
corresponding  to  different  effective  field  value. 

CONCLUSION 

In  sum,  a  new  MOSFET  gate  injection  current 
calculation  method  which  combines  the  attributes  of  the 

MC  and  ET  lechruques  is  presented.  Calculated  values 
for  gate  current  agree  with  experiment  without  the 
need  of  fitting  parameters.  The  method  allows  for  excel¬ 
lent  resolution  of  the  high-energy  tail  of  the  distribution 
function.  Furthermore,  the  ET-MC  model  requires  less 
CPU  time  than  standard  MC  methods,  and  is  therefore 


appropriate  for  use  in  the  computer  aided  design  of 

semiconductor  devices,  especially  after  establishing  a 

database  of  distribution  functioru. 
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Figure  4  Duinbulion  functions,  calculated  using  ET-MC 
iiKtliod,  at  channel  coordinates  I  (>45piii  anu 

I  061prr.  which  correspond  to  elTcc  vc  field  values  of 
202  KV/cm  and  237  K  Wcni  respectively. 


Figure  5  Values  calculated  for  gate  currcni  using  the  presented 
ET-MC  nKthod  and  Richardson’s  eijuaiion  arc 
compared  vith  experiment  The  asterisks  arc 
calci  laied  using  ET-MC  method;  the  open  ciales  and 
'x'  stgns  arc  calculated  with  Richardson's  ei|uation; 
the  solid  lines  are  cxpenmental  values 


